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Abstract 

We formulate a kinetic theory of self-interacting meson fields with an aim to describe 
the freezeout stage of the space-time evolution of matter in ultrarelativistic nuclear 
collisions. Kinetic equations are obtained from the Heisenberg equation of motion for 
a single component real scalar quantum field taking the mean field approximation 
for the non-linear interaction. The mesonic mean field obeys the classical non- linear 
Klein-Gordon equation with a modification due to the coupling to mesonic quasi- 
^sO \ particle excitations which are expressed in terms of the Wigner functions of the 

quantum fluctuations of the meson field, namely the statistical average of the bilin- 
£f~^ ■ ear forms of the meson creation and annihilation operators. In the long wavelength 

| limit, the equations of motion of the diagonal components of the Wigner functions 

take a form of Vlasov equation with a particle source and sink which arises due to 
the non-vanishing off-diagonal components of the Wigner function expressing co- 
herent pair-creation and pair-annihilation process in the presence of non-uniform 
condensate. We show that in the static homogeneous system, these kinetic equa- 
^ tions reduce to the well-known gap equation in the Hartree approximation, and 

hence they may be considered as a generalization of the Hartree approximation 
method to non-equilibrium systems. As an application of these kinetic equations, 
we compute the dispersion relations of the collective mesonic excitations in the 
system near equilibrium. 



Key words: meson condensate; Vlasov equation 



1 Introduction 



Theoretical study of the space-time evolution of matter in high energy nu- 
clear collisions has a long history since the early pioneering works of Fermi 
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and Landau employing thermodynamics and fluid mechanics More recent 
works [2| have been motivated by the prospect of studying a new form of 
matter experimentally by means of very energetic collisions of heavy nuclei as 
those currently underway at Brookhaven with Relativistic Heavy-Ion Collider 
(RHIC) and planned at CERN with Large Hadron Collider (LHC). The data 
taken from the RHIC experiments has shown existence of a strong anisotropic 
(elliptic) collective flow of hadrons in non-central collisions (3L indicating early 
local equilibration of dense matter produced by the collision |4|. This supports 
the hydrodynamic picture of matter evolution and gives us a hope to learn 
information of the equation of state of dense matter from systematic analyses 
of the data. 



Much attention has been paid to the early stage of matter evolution where 
a dense plasma of unconfined quarks and gluons has been expected to be 
formed through some complex non-equilibrium processes [BJ . One anticipates 
also a breakdown of the hydrodynamic behavior of the system when it is 
diluted sufficiently and the collision time exceeds the characteristic time scale 
of expansion. This stage of the matter evolution is usually refered to as the 
freezeout stage. The aim of this work is to present a kinetic theory which is 
designed to describe the freezeout stage of the expanding hadron gas. 

The freezeout of the expanding hadronic matter would proceed in several 
steps; chemical freezeout of the relative abandunce of hadron spieces may 
occur before the kinetic freezeout of the momentum distribution of hadrons. 
Observed relative abandunce of hadrons fit very nicely to a simple picture that 
the chemical freezeout occurs at a certain temperature and baryon chemical 
potential @,[3|. Some of the other observables, such as two particle momentum 
correlations, the pion analogue of the Hanbury Brown and Twiss two photon 
intensity interferometry, may be very sensitive to the dynamics of the final 
kinetic freezeout as indicated by the "HBT puzzle" found in the recent data 
analysis Q. 

If the initial state of the dense matter formed in the collision is a plasma of 
deconfined quarks and gluons, freezeout of the color degrees of freedom should 
proceed these processes when the quark-gluon plasma hadronizes. Chiral sym- 
mmetry an approximate global symmetry of QCD which becomes exact 
in the limit of vanishing quark masses and is considered to be broken sponta- 
neously in the QCD vacuum, may also play an important role in the freezeout 
dynamics. 

As a dense hadronic matter formed by ultrarelativistic nuclear collision is di- 
luted by the expansion, one expects that the system undergoes a phase transi- 
tion associated with the spontaneous breakdown of the chiral symmetry which 
is restored temporalily after the collision by the formation of a quark-gluon 
plasma. As the quark-gluon plasma hadronizes and the system turns into the 
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confining phase, the system would gradually develop a vacuum chiral conden- 
sate and remaining excitaions would expand in the influence of the growing 
chiral condensate. This physical picture has been elaborated in terms of a 
classical equation of motion for the chiral condensate; effects of excitations 
were described in terms of statistical fluctuations in the classical fields 
The fluctuating condensate described by classical pion field has been termed 
Disoriented Chiral Condensate (DCC), emphasizing the symmetry aspect of 
the problem [10]. What was missing in these classical treatments of the meson 
fields is the existence of particle excitations in addition to the condensate. In- 
clusion of particle excitations requires the quantization of the fields. The effect 
of the quantum fluctuations of meson fields has been studied by Tsue, Vau- 
therin and one of t he p resent authors (TM) [ll| by the functional Schrodinger 
picture formalism [l2|, [l3| . 



The physical picture we have described is very similar to what happens when 
the dilute gas of magnetically trapped atoms of alkali metals cools down by 



evapolation [14j . Some of the atoms condense into the lowest single particle 
level in the trapping external potential forming a Bose-Einstein condensate. 
The dynamics of such a system may be described by the coupled equations 



of motion for the condensate, or the Gross-Pitaevskii equation [15] , and the 



kinetic equation, or the Boltzmann-Vlasov equations, for the phas e sp ace dis- 



tribution of the excitations in the presence of the condensate [16|, [17 ■ 



In this paper we will show that a similar set of equations can be derived for a 
system of interacting mesons described by the relativistic quantum field theory 
by the mean field approximation. This approximation corresponds to a neglect 



of all correlations in the system [la ]. We will make no attempt in this work to 



justify this approximation and leave it as an open problem for further study. 

In this context we note that field theoretical derivation of a Boltzmann-type 
kinetic equation has been given by many others 19j. Most of these works fo- 
cuses on the derivation of the collision terms. The present work is distinguished 
from these works in the emphasis of the role of the mean field in the evolution 
of the system coupled with the quasi-particle excitations in the same spirit as 
in 13|, yjj . We omit the effect of the quasi-particle collisions in this work. This 
procedure may be reasonable, at least as a first step, for describing dynamical 
aspect of the freeze-out process: even in the absence of collisions, interactions 
between the quasi-particles and the evolving condensate would affect the fi- 
nal particle distribution. But the present approach is not adequate, however, 
for the early thermalization problem, where the collision terms play essential 
role 0] . We adopt the standard Heisenberg picture instead of the functional 



1 We note that the mean fields may also play important role in the thermalization 
problem, for example via the non-Abelian analogue of the plasma instability [20f| . or 



by the anomalous enhancement of collision rates in a turbulent plasma evolution 2l|] . 



3 



Schrodinger picture used in [13|, since we found it more straightforward to 
see the connection to familiar semi-classical kinetic equation. In this first of a 
series of papers we shall use one-component real scalar field model interacting 
via a 4 self-interaction term in order to concentrate on the presentation of 
the basic features of the theory. Analysis with mult i- component scalar model 
with O(N) symmetry will be deferred for the forthcoming paper [22J. 

In the next section we formulate a quantum kinetic theory for quantum scalar 
meson field, starting from the Heisenberg equations for quantum field of one- 
component real scalar. The mean field approximation replaces the products of 
quantum fields by the products of classical fields and the statistical average of 
bilinear forms of the quantum fluctuations. The latter is expressed in terms of 
the Wigner functions which is reduced to the single particle distribution func- 
tion in the classical Boltzmann equation. Our theory contains another forms 
of the Wigner functions which have no classical counter parts and arises due to 
the coherent pair creation and annihilation processes in non-uniform systems. 
Only in uniform systems, these off-diagonal components of the fluctuation can 
be eliminated by suitable redefinition of the particle mass. Appearance of the 
off-diagonal Wigner functions is reminiscent of the anomalous propagators in 



the microscopic theory of superconductivity [231. \24. I25lj . Similar structure also 
appears in the theory of Bose-Einstein Condensate 17J. Some details of the 
mean field calculation is given in Appendix. 

In section 3, we apply our method to uniform systems and show that each 
mode characterized by the particle momentum obeys non-linear forced oscil- 
latory motion. In section 4 we show that in equilibrium these kinetic equations 



are reduced to a gap equation [26[ which determine the equilibrium amplitude 
of the condensate and the mass parameter as a function of the temperature. In 
this paper, we ignore the effect of the divergent vacuum polarizations which re- 



quires a subtle renormalization procedure in the mean field approximation [27 



The solution exhibits characteristic features of the first order phase transition. 

In section 5, we study slowly varying non-uniform systems. Taking long wave- 
length approximation, the equations of motion of the diagonal components of 
the Wigner functions are reduced to a Vlasov equation in a form generalized 



by Landau for quasi-particles excitations in quantum Fermi liquid [28|,|29( with 
the quasi-particle energy given in the mean field approximation. Non- vanishing 
off-diagonal components of the Wigner functions generate extra terms in the 
Vlasov equation which may be interpreted as particle source and sink terms. 

In section 6, we compute the dispersion relations of excitations of the sys- 
tem near equilibrium. Solving the coupled kinetic equations by linearizing 
the equations with respect to small oscillatory deviations from equilibrium 
solution we obtain dispersion relation of the excitation modes in the system 
near equilibrium. We find the continuum of quasi-particle excitations in the 
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entire space-like energy-momentum region in addition to the continuum in 
time-like region due to the (thermally induced) pair creation. We found that 
in the low temperature phase the meson pole shifts due to the coupling to the 
quasi-particle continua. The effective meson mass vanishes at the edge of the 
spinodal instability line of the first order transition^] We also found that the 
coupling of the off-diagonal components of the Wigner function plays an im- 
portant role to prevent appearance of undamped tachyonic sound mode which 
propagates with a velocity greater than that of light. 

A short summary of the paper is given in section 7 with remarks on the 
remaining problems. 



2 Kinetic equations for the meson condensate and quasi-particle 
excitations 



In this section we derive quantum kinetic equations which describes the time 
evolution of the meson condensate coupled with mesonic quasi-particle exci- 
tations. We use the natural unit h = c = 1 throughout this paper. 



2. 1 Quantized real scalar field in the Heisenberg representation 



We first take a simple model of a self-interacting real scalar field in the Heisen- 
berg picture. The Hamiltonian is given by 



H 



dr 



V + i(V0) 2 + V[0] 



where 



(2) 



Here the scalar field and its canonical conjugate momentum field 7r are 
quantized by the equal-time commutation relations: 



[0(r,t),7r(r',t)] =i5(r - r) 
0(r,t),0(r / ,t)l=[vr(r, t),n(r' ,t)] = 



(3) 
(4) 



2 Softening of the a meson mode has been studied first by Hatsuda and Kunihiro 
as a precursory phenomenon of the chiral phase transition using the Nambu-Jona- 
Lasinio model with quark fields. [3ol| 
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The Heisenberg equation of motion of the quantum field 0(r, t) is given by 
^ = -i[lH]=n(v,t) (5) 



while the equation of motion of the canonical conjugate field 7r(r, t) becomes 
^ = H] = (V 2 - m 2 )0(r, t) - ^A0 3 (r, t) (6) 



Elimination of the field momentum n from these equations yields a modified 
Klein-Gordon equation for the quantum scalar field <j>: 

□0(r, t) + m 2 0(r, t) = "^A0 3 (r, t) (7) 



where □ = d 2 /dt 2 - V 2 . 

2.2 Density matrix and Gaussian Ansatz 



We are interested in the time evolution of the system described by the density 
operator 

p = EI*->p.<*-I ( 8 ) 



where {|^ s )} are a set of normalized wave functions and p s is the probability 
distribution for a mixed state described by this density matrix so that it 
satisfies 

T,Ps = 1 (9) 



This density matrix can be also expressed in terms of some complete set \a) 
of the wave functions of our Hilbert-Fock space as 

P = £|a)M/3| (10) 
where 

Pa p = J2(a\y s )Ps(y s \P) (11) 
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All physical information of the system to be described are contained in a 
specific form of the density matrix. In thermodynamic equilibrium, the density 
matrix is given by 

p cq = (12) 

with 



Z = tr 



(13) 



where F(T) gives the Helmholtz free energy of the system at temperature 
T = 

In the Heisenberg picture the density matrix is time-independent since the 
wave functions are time- independent; all time dependence arises from the time 
dependence of an operator: 

(6(t)) = tr \6{t)p] (14) 



For example, we define the classical condensate fields by the statistical average 
of the quantum fields, 



c (r,t) = (0(r,t)>, (15) 
7r c (r,i) = <7r(r,t)) (16) 

In the following we take a Gaussian Ansatz for the density matrix: 

(0 n (r, £)) = for odd integer; n (17) 

and 

~ n} ~ 

(0 n (r,t)) = — — (0 2 (r,t)) m for even integer; n = 2m (18) 

where 0(r,t) and 7f(r, t) are shifted field operators defined by 

0(r,t) = 0(r,t)-0 c (r,t) (19) 
7f(r,*) = 7r(r,t)-7r c (r,i) (20) 

The shifted field operators obey the same equal-time commutation relations 
as the original fields: 
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(r,t),7r(r',t)] =i5(r - r) 



(21) 
(22) 



We will show that this choice of the density matrix is a non-equilibrium gener- 
alization of the Hartree approximation in equilibrium. It has been alternatively 



introduced in the functional Schrodinger representation [13|, [11 



2. 3 The Wigner functions 



We introduce the one-particle Wigner function by 
F(p,k,t) = (a P+k/2 (t)a p _ k/2 (t)) 



(23) 



Here the particle creation and annihilation operators may be defined in terms 
of the Fourier transforms of the shifted fields 



as 



(24) 



ap(t) 



u p (f)p(t) - i7f_ p (t) 



with 



= \JP 2 + ^ 2 



(25) 
(26) 

(27) 



These relations are rewritten as: 



<KM) = ]>>^ r -— k(t) + ol p (t) 

P V Zuj p 



The quantization rules (121] 1221) are transcribed to 



(28) 
(29) 



K(t),al,(t)]= 8 PtP > 
K(t),a pl (t)} = [aUt),al(t)} = 0. 



(30) 
(31) 
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with which we may interpret a p (at,) as annihilation (creation) operator of 
"particle excitation" with momentum p. 

We also introduce the following other forms of the Wigner functions: 



G(p, k, t) = (a_ p _ k / 2 (t)ap_ k / 2 (t)), (32) 

G(p,k,i) = (a p+k/2 (i)al p+k/2 (i)>, (33) 

F(p, k, t) = (a_ p _ k/2 (t)a t _ p+k/2 (t)), (34) 

These Wigner functions are not independent but are related to each other. 
The complex conjugate of the Wigner functions are given by 

F*(p,k,i)=F(p,-k,f), (35) 

F*(p,M)=F(p,-M) (36) 

G(p,k,t)=G*(p,-k,t) (37) 



where the asterisk (*) stands for the complex conjugate. The commutation 
relations imply also that G(p, k, t) and G(p, k, t) are even functions of p 



G(p,M) = G(-p,M) 
G(p,k,t) = G(-p,k,t) 



(38) 
(39) 



and 



F(p,k,t) =F(-p,M) + &,o 



(40) 



The four Wigner functions may be grouped together to form a matrix form of 
the Wigner function 

, F(p,k,t) G(p,k,t) , 
W(p,k,t)=\ ^ _ lP ' (41) 
G(p,k,t) F(p,k,t) 



The appearance of the "off-diagonal" components of the Wigner functions is 
reminiscent of the anomalous propagators in the BCS theory of superconduc- 



tivity which arises due to the presence of fermion pair condensate |23l . |24| . |25 



Note that our definition of the matrix components of the Wigner function is 
slightly different from those for the propagators. 

We write the Fourier transforms of the Wigner functions as 
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/(p, r, t) = J2 e" ik - r F(p, k, t), /(p, r, t) = £ e- lkr F(p, k, t), (42) 

k k 

<?(p, r, t) = £ e- lk ' r G(p, k, f) , 5(p, r, t) = E e~* r G(p, k, t) (43) 

k k 

The ( 135]) and (|36|) imply that /(p, r, £) and /(p, r, t) are real functions and 
are related to each other by 

/(p,r,t) = /(-p,r,t) + l, (44) 



while ( 1371) implies that g(p,r,t) and g(p, r, t) are complex conjugate to each 
other: 



0*(p,r,t) = S(p.r,t) (45) 

Here we have chosen the particle "mass" \x to be different from the mass pa- 
rameter m in the original Hamiltonian. Physical particle mass for interacting 
fields is generally different from the mass parameter in the Hamiltonian or 
Lagrangian due to the effect of interaction, e.g. renormalization with or with- 
out the spontaneous symmetry breaking. It may also depend on the physical 
conditions described by the statistical average with the density matrix p. Since 
we are interested in non-equilibrium time-evolution of the system where the 
physical particle mass may not have a definite meaning in the intermediate 
states, we consider here p just as a parameter to be chosen at our discretion, 
for an appropriate choice of the initial conditions specified by the Gaussian 
density matrix. A different choice of this mass parameter would give different 
definition for "particle excitations" ; the Wigner functions thus depend on the 
particular choice of the mass parameter. 

Suppose we take a different particle mass p! to define the particle creation and 
annihilation operators 



b p (t) = - 1 L= L p > p (t) + in _ p (t)l (46) 
bl(t) = ~j^= kVp(t) - Mf-p(*)l (47) 



2< 



with 



These new particle creation and annihilation operators should also obey the 
commutation relations, 
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mo, =w 

[b v (t)Mt)] = [bUt),iL(t)] = o. 



(49) 
(50) 



so that they are related to the original ones by the Bogoliubov transformation: 



bp = cosh a p a p + sinh a p aL p 



bL p = sinh a p a p + cosh 



Gip CI _ p 



(51) 
(52) 



where the real parameter a p is determined by requiring that they describe the 
same fields: 



0M) = E 



D ipr 



i p (t) + al p (t)} = £ 



b p (t) + bl p (t) 



(53) 



and this gives 



u' 



UJ r . 



(54) 



or 



a. 



2 l0g \ u„ , 



1 + A 

4 ° g Vp 2 + ^ 2 , 



(55) 



The new Wigner functions defined by replacing the creation and annihilation 
by the new ones are related to the original Wigner functions by: 



W'(p,k,t) 



V( P ,k,t) g"(p,M 

v G'( P ,k,t) F'(p,M 
M(p + k/2) W(p, k, t) M (p - k/2) 



(56) 



where 



M(p) 



coshap sinh q; p 
sinhcip cosha p 



(57) 
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with 

n = \ \, t 2 = I i . , I ' | m 



For a small change of the mass parameter fi — > fi + 5/u the Wigner function 
will change by 






5W(p, k, £) = riW(p, k, t)<kt p+k /2 + W(p, k, t)ri(5a p _ k / 2 
U(p + k/2)2 + ^) + 4 ((p - k/2)2 + fj?) , 

In particular, 



fibfi 

(59) 



5F(P ' k ' t} = (iOpTWT^) + 4(( P -k/2) 2 + ^) J * (60) 



In uniform equilibrium system the mass parameter may be chosen to "diag- 
onalize" the one-body mean field Hamiltonian. As will be shown later, this 
procedure will lead to the well-known gap equation, a self-consistency condi- 
tion to determine fx. If the system is slowly changing in time, one may still 
use such procedures adjusted to slowly varying quasi-equiliblium conditions, 
introducing a time dependent effective mass as a dynamical parameter to de- 
scribe such adiabatic process. For such calculations, the relation (160|) may be 
used to describe the adiabatic change of the Wigner functions by the change 
of the mass parameter. 

In more general non-equilibrium situations as we expect to encounter at the 
freeze-out stage of expanding matter, however, there may be no such appro- 
priate condition to determine the mass parameter. The situation could even 
be worse: if the system goes through unstable state with respect to small fluc- 
tuation of the field, then the adiabatically determined mass parameter would 
become pure imaginary, reflecting the extremum, instead of the minimum, of 
the effective potential. In such case we may keep the value of fi at a certain 
real value reflecting initial conditions. The instability would then show up 
as appearance of a growing solution to our kinetic equations; which would 
eventually be stabilized by the non-linear interaction. 

To extract the physical information, such as the particle distribution in the 
final asymptotic state, we should use the Wigner function defined with the 
physical particle mass in the vacuum. But these asymptotic physical Wigner 
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functions may be calculated from the Wigner functions with a different choice 
of the mass parameter by the relation (|56|) . 



2.4 Equation of motion in the mean field approximation 



The equation of motion of the classical mean field c (r, t) is obtained by taking 
the quantum statistical average of the field equation © . With the Gaussian 
Ansatz for the density matrix, we find 



□0 c (r,t) + m 2 c (r,t) 



1 

3! 



A 



^(r,t) + 3(0 2 (r,t))0 c (r,t) 



(61) 



This equation corresponds to the non-linear Schodinger equation (also called 
the Gross-Pitaevskii equation 14] ) in the theory of Bose-Einstein conden- 
sates. So we may call this equation non-linear Klein-Gordon equation. The 
non-linearity arises due to the self-interaction of the classical field </> c (r, t) 
(condensate) and also due to the interaction with fluctuations (</> 2 (r, t)) which 
also depends on <p c (r,t) implicitly. The latter may be interpreted as due to 
"particle excitations", since the fluctuation can be expressed by the Wigner 
functions as 



e i(-p+p')-r 

' t )) = E ^~ /rj— (( a -p + a P )(V + a -p')) 

P ,p' V Zu pV Zuj p' 



p j(-p+p')-r 

e4 



p,p 



p + p' 



F{ — - — ,p-p,t) +F( — - — ,p-p,t) 



P - P ',t) + G(P±P-, P - P ',t) 



(62) 



The time-evolution of the classical mean field 4> c (r, t) is thus coupled with the 
time-evolution of the Wigner functions. 

To derive the equation of motion of the Wigner functions, we need to compute 
the time-derivative of the bilinear forms of the operators a p (t) and a p (t) which 
in turn requires computation of the commutators of these operators with the 
hamiltonian. We decompose the original hamiltonian as 



H = Hq + H\ + H2 + H3 + H4 



(63) 



where Hq is the classical hamiltonian obtained from H by replacing the quan- 
tum fields by their classical expectation values and Hi contain the z-thrth 
power of the quantum fluctuation (or a p and a P ). A straightforward calcu- 
lation yields 
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#1 

H 2 



dr 
dr 



7T C 7T + V0 c V0 + m 2 (f) c 



A 
3!' 



^ + ^(V0) 2 + in 2 2 M) + ^> 2 M)0 2 (r,t) 



A 



H 3 = - J dr^r,^^^) 
H A = ± j dr^(r, t) 



(64) 
(65) 

(66) 
(67) 



The commutators of bilinear forms of a p (t) and a p (t) with if x vanish and the 
commutators with H 3 would give either a linear term or the third power of 
the fluctuation, both of which may vanish when taking the average with the 
Gaussian density matrix. What remain to be computed are then the commu- 
tators with if 2 and with H4. They will give either the bilinear form of a p (t) 
and a p or the fourth power of the fluctuations. The Gaussian average of the 
resultant equations of motion of the bilinear field operators would give the de- 
sired equations of motion of the Wigner functions. Details of this computation 
is given in Appendix A. 

The resultant equation of motion of the Wigner functions may be obtained 
more easily by introducing the mean field Hamiltonian defined by 



mf 



= /dr 
= /dr 



V + i(V0) 2 + im 2 2 + in(r,t)0 2 
^ 2 + i(V0) 2 + ^ 2 2 + ^An(r,t)0 2 



(68) 



where 



n(r,t) = ^(0 2 (r,t) + (0 2 (r ) t))) 



(69) 



and 



AII(r,i) = II(r,t) +m 2 - /J 



(70) 



In the momentum representation this mean field Hamiltonian may be written 

as 

ii _\^,,„t„ 1 1 VAn (a p + a 1 lp)(a-p-q + a P+q ) , , 

p ^ P,q VHVH+q 
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where 

AII q (t) = J rfre- iq ' r (n(r, t) + m 2 - /i 2 ) = II q (i) + (m 2 - /z 2 )^ (72) 



The commutator of a bilinear operator product of a p and ctp with this mean- 
field Hamiltonian is given by 

[<a P2 , # mf ] = -Ki - ^Jat iQp2 - £ An q • (a - p / ^ a ^ P2 

q V Za; Pi+qV Ztt; P2 

+EAnq .<Kp±|kj 2 ) (73) 

We show in Appendix that the quantum statistical average of this commutator 
with the Gaussian density matrix gives precisely the same result for the same 
statistical average of the commutator with the original Hamiltonian: 

([a Pl a P2 ,# mf ]> = ([a pl a P2 ,H]) (74) 



Therefore one can compute the equations of motion of the Wigner functions 
using this effective Hamiltonian, 

d 

^ F (P> k > *) = ([ a P +k/2 a P-k/2, H}) = ([a P+k/2 a p _ k/2 , H mi ]) (75) 

Using this we find, 
d 

2 ^ F (P' k > *) = -(^P+k/2 - u p -k/2)F(p, k, t) 

F(p + q/2, k + q, t) + G(p + q/2, k + q, t) 



E^Hq 

q 

+ E An q 



/ 2cJ p+ k/2 \J 2dJ p+ k/2+q 

F(p - q/2, k + q, t) + G(p - q/2, k + q, t) 

^ 2o;p_k/2 \J 2^p-k/2-q 

(76) 



Equations of motion of other three Wigner functions G(p,k, t), G(p, k, t), 
F(p, k, t) can be also computed from the commution relations of the prod- 
uct operators a Pl a P2 , a pi a p2 , a Pl a p2 , with the mean field Hamiltonian H m f, 
respectively: We obtain 



d 

i— G(p, k, t) = H+k/2 + w p _ k / 2 )G(p, k, t) 
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+ ^ An G(p + q/2, k + q, t) + Fjp + q/2, k + q,t) 



+ £An 



2co> p+k /2 y 2^p+k/2+q 

G(p - q/2, k + q, t) + F(p - q/2, k + q, t) 



q 



q y2u;p_ k /2y 2u; p _ k / 2 -q 

(77) 

d - 

^C(p, k, t) = -(u; p+k/2 + Wp_ k / 2 )G(p, k, i) 

_^ An g(p + q/2, k + q, t) + F(p + q/2, k + q, t) 



q y 2l^p+ k /2y 2^ p+k /2+q 

_^ A n ^(P - qA k + q, t) + F(p - q/2, k + q, t) 



2co>p_ k / 2 Y 2co , p _ k /2-q 



(78) 



and the equation of motion of F(p, k, t) can be obtained from (1751) by the 
substitution p — > — p: 



- 

*t^(p, k, t) = (^ p+k /2 - Wp_ k / 2 )F(p, k, t) 

i ^An q f (p + qA k + q,i) + G(p + qAk + q,t) 

2^ p+k /2 Y 2^ p+k /2+q 



^ An F(p - q/2, k + q, t) + G(p - q/2, k + q, t) 



2o; p _ k /2 Y 207p_ k /2-q 

(79) 

These equations form a closed system of coupled differential equations with 
the non-linear Klein-Gordon equation (16TT) which may be rewritten as 

□0 c (r,t) + (/i 2 + An(r,t))0 c (r,t) = -^(r,t). (80) 

These four equations of motion of the Wigner function may be combined into 
a single matrix form as 



d 

i-Q^W(p, k, t) = -uj p+k/2 r 3 W(p, k, t) + w P -k/2^(p, k, i)r 3 



r 3 (l + r 1 )^(p + q/2,k + q,t) 



2^p+k/2 Y 2Wp+ k /2+q 



+ ^ A n .^(p-q/ 2 ' k + ( i^)( 1 + r i) 7 "3 



2w p -k/2 Y 2u;p_ k /2-q 
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(81) 



3 Uniform system 



For a uniform system, we expect that the classical mean field and the self- 
energy become functions only of time: 

c (r,t) =Mt), n(r,t)=n (t) (82) 



Thus the non-linear Klein-Gordon equation (|80|) becomes 

0o(r, t) + (fi 2 + AU (t))Mt) = (83) 



and the mean-field Hamiltonian is reduced to 

H mi = Y,^P a W + AIIo(t) Y,( a p + a - P )( a -P + Op) ( 84 ) 



In this case the Wigner functions contain non-vanishing components only for 
the diagonal elements (k = 0) so that they may be written as 



F(p,k,t) = F (p,t)(J(k) (85) 

G(p,M) = G (p,t)<S(k) (86) 

G(p,k,t) = G (p,t)5(k) (87) 

F(p,k,t) = (F (-p,t) + l)5 k ,o (88) 

Then the equations of motion of the Wigner functions become 



i j^b(p, t) = ^ (G (p, t) - G (P, 0) (89) 

i j^o(p, f) = ^ (G (p, t) - G (p, t)) (90) 
i^G (p, t) = + G (p, t) 

All °(Fo(p,t) + F (p,t)) (91) 



2u; p 

i^G (p, t) = - (^ p + ^) G (p, t) 
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Ail " 'F (p,t)+F (p,t)) (92) 



2u p 



These coupled equations have a time-independent solution of the form 

Go(p) = <5 (p) = c (Fo(p) + F (p)) (93) 

with 

An n ( An \ _1 , x 

2lu p + °-\ (94) 



UJ p \ UJ p 



Thus if we define the particle mass /i so as to satisfy Alio — 0, then two 
components of the Wigner functions Go(p) and Go(p) vanish for all p. In 
this case, F (p) may be interpreted as the momentum distribution of physical 
particle excitations. 

To find time-dependent solutions, we write 



F ± (p.t)=F (p,t)±F Q (p,t) (95) 
G ± (p.t)=G (p,t)±G (p,t) (96) 

and rewrite the equations as 



G_(p,t) 



• 9 TP ( A An 0(^) 

<|F_(p,t) = 
d„ , \ ( AIIo(t)' 

,-G_(p,t)= 2. p + — A2 



G_(p,t) 



An 



G + (p,t) + -F + (p,f) 



(97) 
(98) 
(99) 

(100) 



These equations are to be solved together with the non-linear Klein-Gordon 
equation fl83l) with 



An (t) = ^ 



^) + E F+(p,t l +G+(p,t) +^ 2 -M 2 



2uv 



(101) 



From these equations we can derive second order differential equation: 
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d_ 2 
dt 



q2 



/All N 2 



o ^ An / An \ . 

F+(p, t) - —— 2u p + G+(p, t) 



+i 



1 dAn (t) 
(it 



GL(p,t) 



(102) 



An n / An 

2cj p + 



^'F + (p,t)- (2^ P + ^) G +(P ,t) 



UJr 



. 1 dAIIo(t) 

-z — G-{p,t) 

uj p at 



(103) 



with 



dAIIo(t) 



fro(*)0o(*) + £ G -(P-*) 



(104) 



We observe that each of these coupled equations looks like an equation of a 
forced oscillator. For example the frequency of the oscillator F + is given by 
■^2^ while the frequency of the oscillator G + is given by 2co> p + so that 
G + oscillates much rapidly than F + . The coupling of these oscillators may 
produce an interesting interesting effects which may be studied by numerical 
integration. 



4 Statistical equilibrium 



In statistical equilibrium, we have the Bose distribution 



^o(p) = / cq .(p) = ^3 T (105) 
together with 

Go(p) = 2c q .(p)=0. (106) 



where (3 = 1/ksT is the inverse temperature. The condition of vanishing An 
implies from f!72|) 



~ (ft + (0V) + m 2 - /i 2 = 0. (107) 
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where the thermal fluctuation of the quantum field is given by the relation 
(1521) as 



2 )eq. = /eq.(P) 



p 



(108) 



On the other hand, the condensate amplitude 0o i n equilibrium should also 
satisfy the static non-linear Klein-Gordon equation, 



^ + 3(0 2 ) eq .0 o 



(109) 



From these two conditions, we find, for m 2 < 0, 



-2m 2 -X^— /eq.(P 



p ^p 



(110) 



and 



A0Q = 3/i 2 



;ni) 



We note that, since the thermal distribution / eq . depends on the mass (gap) 
parameter \x through u> p = ^/p 2 + /i 2 , the equation (11 10p determines the mass 
parameter \x self-consistently as a function of the temperature. This equation 
is called the gap equation 26|, |27| • As the temperature increases the mass gap \x 
and the condensate amplitude 0o decreases and vanishes at high temperature. 



To find the behavior of the mass parameter as a function of the temperature we 
can use the following formula 26|, |34| for the integral of the bose distribution 
function: 



J2— Zeq-(P) 



P ^P 



dp 



(2vr) 3 J u p e^f 3 - 1 2vr 2 
where the dimensionless function I_ (x) is given as 



(112) 



r(2), 



X = 



ifc 2 difc 



1 



7T 7T 1 2 X 
.T X In 

6 2 4 4tt 



^ + <?(* 6 ) 

OZ7T 



(113) 
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(7 = 0.57721 • • • is Euler's number) . The gap equation (lllOp has a solution 
H = at T c determined by 

24 

(k B T c ) 2 = -—m 2 (114) 



However the solution exhibits the behavior of the first order transition due 
to the non-analytic behavior of the function I^\x) This is a generic 

feature of the mean field approximation [27| which may be a theoretical artifact 
and may not survive the inclusion of correlations missing in the mean field 
approximation. An improvement of the mean field approximation has been 
proposed in [35j and later it has been shown to lead to the second order 
transition in two-loop approximations 36 . 
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Fig. 1. Temperature dependence of the mass parameter \x for A = 1 (dashed line) 
and A = 10 (solid line) 

In Fig. 1 we plot O as determined by (11091) as a function of temperature. It 
shows a behavior characteristic of the first order phase transition: there is a 
region T c < T < T where there are three solutions of (11091) , one at 0o — and 
other two at = 0i ^ as well as = O 7^ 0. One expects that two solutions 
= 0o an d = 0o correspond to two local minima of the effective potential 
V(4>) while the other one = 0i corresponds to the local maximum of V (0). 
Indeed we can show this explicitly by constructing the effective potential V(4>) 
by demanding that the local extremal condition dV((p)/d(f) = coincides with 
the condition (1X091) . 

V(J>) = Jd<P (m 2 + ^0 3 + ^E (H5) 
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The integration of the first two terms are trivial. Beside irrelevant integration 
constant, which we may choose to be zero, we obtain 

V{</>) = lm 2 2 + A 4 + A J Lgti (n6) 



To carry out the remaining integral over <\> we demand that the mass parameter 
H depends on <ft through the relation 

V? = m> + \ (0 2 + (0 2 ) cq .) = m 2 + \ (> + £ ^) (117) 



which appears to be equivalent to the condition (11071) . but we assume that 
this relation holds not only for the equilibrium value of (j), namely at <p = <fio, 
but also for any value of <fi. Noting 




which results from (11171) . we transform the integration variable from to /i 2 : 



A f J^^ST^ /Cq-(P) _ 1 f d[ /(l X 9 y /eq.(p) A /eq.(p' 



2 



p' 

2 



2 y M v ^ P 8 V p ^ / 



+ const. 

(119) 

Remaining integral over fi 2 can be carried out as 



2 /eq.(P) 1 m / a. . M V- P 2 



- £ In [1 - «p(-/H)] = - E f-/e q .( P ) (120) 



where in deriving the final expression we performed integration by part in p 
and omitted again the irrelevant integration constant. 

To sum up, we find the following expression for the effective potential: 

VW - imV + ^ - E £/.-(p) - I (E + const. (121) 



where <p dependence of the last two terms are implicitly given from the 0- 
dependence of \i. We note that this result coincides, in the neglect of the 
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renormalization effects due to the meson mass shift, with the effective potential 



obtained by Amelino-Camelia and Pi [37| using the Cornwall- Jackiw-Tombolis 



(CJT) composite operator effective potential formalism [38], We note that our 
mass parameter fi corresponds to the variational mass parameter M of the 
CJT potential. The third term may be interpreted as the pressure of an ideal 
gas of the quasi-particles obeying the dispersion e = u p = y/p 2 + /i 2 . 

The effective potential given by the formula (11211) takes in general a complex 
value since fi 2 becomes negative at low temperature and at small value of <fi. At 
low temperatures the potential develops a second local minimum at non-zero 
value of <fr, which becomes the absolute minimum below certain temperature 
Ti (T c < T\ < T ) where the first order transition takes place in equilibrium and 
the order parameter O jumps discontinuously from to 0o as the temperature 
is lowered. 



5 Slowly varying system: the Vlasov equation 



In the presence of inhomogeneity the Wigner functions acquire non-vanishing 
elements with k ^ 0. If we assume that this inhomogeneity are due to the 
long-wavelength fluctuations in the system, 

k,q<^p (122) 



In this case we can obtain a familiar form of the Vlasov equation [39( from the 
equation of motion of the Wigner function by the procedure usually referred 
to as the gradient expansion 4^ . 



To show this we make the following approximations: 



^p+k/2 + ^p-k/2 - 2u; p (123) 
p • k 

^p+k/2 - ^p-k/2 ^ (124) 



and 



2Wp+k/2 X /2^p + k/2+q 2UJ P V 



if p-(k + q)\ , 

'l- ,. a (125) 



if p • (k + q) \ . 

1 + .a ( 126 ) 



2wp_ k / 2 y / 2^ p _ k /2- q 2uJ P \ u p / 

We also make use of the Taylor expansion of the Wigner functions: 
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F(p ± q/2, k + q, t) ~ F(p, k + q, t) ± ~q • V p F(p, k + q, t) (127) 
G(p + q/2, k + q, t) ~ G(p, k + q, t) + iq • V p G(p, k + q, t) (128) 
G(p - q/2, k + q, t) ~ G(p, k + q, t) - ^q • V p G(p, k + q, t) (129) 



Then the equations of motion (1761) of the Wigner function F(p, k, i) becomes 



|-F(p, k, t) = t^F(p, k,t)+iJ2 • V p F(p, k + q, t) 

Ol LU p q ZLUp 



An q p_(k + q) 

q ^p 

+* E ^ ( G (P- k + q, t) - G(p, k + q, t 

q ZUJ P 



~i E t nqP ' (k ^" q) ( G (P' k + q, t) + G(p, k + q, t 



q 2 ^p W p 



E ■ V P (G(p, k + q, t) + G(p, k + q, t) 



q 4 ^P 



(130) 



On the right hand side, the first term and the second term are disguised forms 
of the drift term and the Vlasov term respectively. 

To obtain more familiar form we make the Fourier transforms (42) and (43) 
of the Wigner functions. Then we find 



|/(p, r, t) = - P- • V r /(p, r, t) + V r (^^-) ■ V p /(p, r, t) 
at uj p \ lujp J 

AII(r,t) p „ „. 

2^ p LUp 2 

.Afl(r,t) / / / 

+t K (g(p,r,t)-g(p,r,t)) 

LLOp 

Afl(r,t) p „ , , 
+ ■ V r (g(p, r, t) + g(p, r, t)) 



2ujp 2 Up 



+V r ^4^^ ) • Vp (g(p, r, t) + g(p, r, t)) 



;i3i) 



Now it is clear that the first term is the drift term which describes the change 
of the particle position by the drift with the velocity v p = p/ou p . The second 
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term can be interpreted as the Vlasov term which represents the change of par- 
ticle momentum due to the continuous acceleration by the velocity dependent 
equivalent potential, 



acting on the particle. The third term appears to be the correction to the 
drift term due to the local change of particle mass //—>// = // + AIl(r, t) 
which causes change in particle velocity. Other three terms are associated 
with the other components of the Wigner function and has no counter parts 
in non-relativistic Vlasov equation. 

Noting that 



(133) 



dp u p 

and 

dU P (r,t) = U p (r,t)p 

Op LUp 2 



the above equation may be rewritten in a more compact form: 

-|/(p, r, t) +V p e(p, r, t) • V r /(p, r, t) - V r e(p, r, t) ■ V p /(p, r, t) = 
iU p (r, t)g-(p, r, t) - ^V p C/ p (r, t) ■ V r g + (p, r, t) 
+-W r U p (r,t)-W p g + (p,r,t) 

(135) 

where 

e(p, r, t) = Up + U p (r, t) (136) 



and 



9±(P, r, t) = g(p, r, t) ± ^(p, r, t). 



(137) 
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The quantity e(p, r, t) plays the same role as the quasi-particle energy which 
appears in the kinetic equation of Landau's Fermi-liquid theory [281. We 
note that because of the relation (145 p . g~{p, r, t) is a pure imaginary function, 
while <7+(p, r, t), /(p, r, t), U p (r,t), and e(p,r,t) are all real functions. 

In the long wavelength approximation, the equation of motion of G(p, k, t) 
becomes 



|rG(p, k, t) = -2iu p G(p, k, t) - i E — G(P, k + q, t) 
at q ujp 



~* E 4^ (^(P, k + q, t) + F(p, k + q, t) 



q 2 ^P 



+? ^ m^ p-^ + qj / F(pj k + q t) _ F(pj k + ^ t) 



An q p ■ (k + q 

q 2 ^P ^ P 

E • V p ( F (P' k + q, t) - F(p, k + q, t 

q 4C ^P 

(138) 



which may be rewritten for the Fourier transforms in a compact form as 



d 

qI9(p, r , t) + 2ze(p, r, t)g(p, r, t) = -iU p (r, t)/+(p, r, £) 

--V p ?7 p (r, t) • V r /_(p, r, t) + X - V r C/ p (r, t) • V p /_(p, r, t) 



(139) 



where /±(p, r, £) = /(p, r, £)±/(p, r, t). Thus the kinetic equation for g(jp, r, t) 
does not look like a Vlasov equation; it takes a form of a simple oscillator 
equation with the oscillator frequency 2e(p, r, t) exposed to the external per- 
turbation created by the particle distribution /(p, r,t). Similar equation is 
derived for g(jp, r, t) which is the complex conjugate of g(p, r, t). 

Taking the complex conjugate of the Eq. (11391) . we obtain the kinetic equation 
for g{p,r,t) = g*(p,r,t): 



d _ 

fr9(P, r > *) - 2 ^(p, r, t)g(p, r, t) = -iU p (r, t)f + (p, r, t) 

-^Vp[/ P (r, t) ■ V r /_(p, r, t) + - V r C/ p (r, £) • V p /_ (p, r, t) 
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Relativistic extension of the general framework of the Landau Fermi-liquid theory 



has been made by Baym and Chin [3l[ and applied by one (TM) of the present 
authors [33] to Walecka's relativistic mean field theory of cold dense nuclear matter 
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(140) 



Adding or subtracting (I139p and (11401) . we find for the real function g+(p, r, t) 
and the pure imaginary function g_(jp,r,t) 



d 

^fl , + (p > r,t) = 22e(p,r,t)p_(p,r,t) - V p C/ p (r,t) ■ V r /_(p,r,t) 
+V r E/- p (r,t).Vp/_(p,r,t) 

(141) 

^g-(P, r, t) = 2ze(p, r, t)g + (p, r, t) - 2iU p (r, t)f + (p, r, t) (142) 



In the long wavelength approximation, we also have 



f(p, r, t) + f(p, r, t) + g(p, r, t) + g(p, r, t) 
p 

= E 77^" [/+(P, r, f) + p+(p, r, t)] (143) 
P Zuj p 

which appears in the non-linear Klein-Gordon equation (IBTl) for the condensate 
and in the self energy term (1691) in the Vlasov equation. 



6 Dispersion relation of the excitations near equilibrium 



We now apply our coupled kinetic equations, consisting of the non-linear Klein- 
Gordon equation and the Vlasov equations, to find a dispersion relation of the 
excitations in the system near equilibrium. The dispersion relations of the 
excitations of the finite temperature system may be computed using the prop- 
agator (two-time) formalism, either with real-times (Schwinger-Keldysh) for- 
malism [42J or with imaginary-times (Matsubara) formalism with subsequent 
analytic continuation k3j. Although such calculations have been performed 



by various authors [351. 144j . |45j numerical results have been presented focusing 
on the time-like region of the spectral functions. We study with our kinetic 
theory the entire range of the (u, k) plane including the space-like excitations 
of the system. As has been shown by one of the present authors 32] the 



We note that [35f] contains a general formula of the meson spectral function cov- 
ering all kinematic ranges including space-like region. Their results are indeed very 
similar to ours, besides that the divergent vacuum loops are included with an elab- 
orate temperature-dependent renormalization procedure and that the "tree-level" 
masses are used in the modified loop calculation, while we ignore the vacuum loops 
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calculation of the zero sound mode in the degenerate Fermi liquid using the 



relativistic Landau kinetic equation[3lj], similar to ours, reproduces the same 



result as the propagator theory in the long wavelength limits. We thus ex- 
pect that our result may also reproduce the long wavelength behaviors of the 
excitations obtained from the propagator theory. 

For this purpose we assume that the distribution functions consist of uniform 
equilibrium term and a small deviation from it: 

/(p,r,t) = / eq .(p) + 5/(p,r,t) (144) 
g(p,r,t) = g cq Xp) + 5g(p,r,t) (145) 

where 



f«i ( p ) = e ^ p /3 _ i and #cq.(p) = 0. (146) 

We also assume the mean meson field consists of uniform equilibrium term 
and a deviation from it 

<f) c (r,t) =0 O + 50M) (147) 



where 0o is determined by the static non-linear Klein-Gordon equation (11091) 
together with the solution of the gap equation (1 11 01) for the mass gap \x. 

Linearization of the non-linear Klein-Gordon equation with respect to 8<p(r, t) 
yields 

□<ty(r, t) + m 2 50(r, t) = — -A [$<ty(r, t) + (4> 2 ) C(l 5<p{v, t) + 6(4> 2 (r, t))<j> 

(148) 

where the thermal fluctuation of the quantum field in equilibrium is given by 

(0V=E-/e q .( P ) (149) 

p Up 



and the deviation 5{4> 2 (r, t)) is expressed in term of the small deviation of the 
distribution function from equilibrium value. 

and the mass of our quasi-particle excitations is determined by the gap equation 
which includes the effect of one-loop diagrams in the language of conventional loop 
expansion. 
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5(4> 2 (r, \Sf(p, r, t) + 5f(p, r, t) + 5g(p, r, t) + <^(p, r, *)' 

P Zuj p 

= E 7T- [ S f(& r ' *) + 5 /(-P' r ' *) + ^(P. r ' *) + 5 9* (P» r > *)] 
P 2cj p 

(150) 

where in deriving the last line we have used 

5f(p,r,t) = 8f(-p,r,t) (151) 
6g(p,T,t)=6g*(p,r,t) (152) 



which follow from (jHj) and (I45j) respectively. 

By the linearization with respect to the small increments, Sf(p, r, t), 5g(p, r, t), 
the Vlasov equations become 



^5f(p, r, t) + V p o; p ■ VJf(p, r, t) - V r 5C/ p (r, *) ■ V p / cq .(p) = 

(153) 

.9 

— (Jflf(p, r, £) + 2zw p 55((p, r, t) = -i26U p (r, *)/«,. (p) 

(154) 

The right hand side of the Vlasov equation (11351) vanishes because AIl(p, r, t) 
vanishes in equilibrium. The fluctuation in the mean field potential SU p (r,t) 
which appears in these equations is related to the fluctuations of the conden- 
sate amplitude and the distribution functions 



6U p (r, t) = 2^"5n(r, t) = ^- (20 o 50(r, t) + 6(4> 2 (r, t)>) (155) 

The linearized Klein-Gordon equation (11481) and the linearized Vlasov equa- 
tions (11531) and (11541) . supplimented by the "constitutive relations" (11551) and 
(11501) . form a closed set of equations to determine a small fluctuation propa- 
gating the system in equilibrium. 

Noting that 5(f>(r, t) and Sf(p, r, t) are real functions, we seek a solution in the 
form: 



6(f)(r, t) = 5<Pe i ^ r ~ UJ+t) + d^e'^-^ , (156) 

5f(p, r, t) = Sf p e^ r -^ + Sfte-^-"^, (157) 

5g(p,r,t) = 5g p e^ r -^ t \ (158) 

Sg(p,r,t) = 5g*(p,r,t) = Sgte^^, (159) 
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where we have introduced the Landau prescription [391 ] . 

lj± = uj ± ie 



(160) 



with a positive infinitesimally small constant e, to set an adiabatic switching- 
on of the fluctuation, namely 5f(p,r,t), 5g(p,r,t) all vanish slowly as t — > 
— oo. The constitutive relations for the fluctuations (11501) and the mean field 
potential (11551) now read 



■r— u)+t) 



,5(0 2 (r,i)> -I- (*/p + */- P + Sg p + fy- p ) e l(k " 
P ZuJ p 



P 2uj p 



(161) 



and 



<5[/ p (r,t) 



2Ur 



2u. 



s <f> + E 7T- ( 6 fp + */-p + s fb + $9- 
P Zu; p 

m* + E «V (*£ + */- P + + s 9 



i(k-r— u>+i) 



P 2 ^p 



-i(k-r— 



(162) 



respectively. 

Using these relations, we find from the linearized Klein-Gordon equation ( 11481) 

2 , l 2 , 2 

-co> + + k + fj, 



^E^-te+^-p+^+^J. ( 163 ) 
^ P ^ P 



and its complex conjugate relation, where we have used the equilibrium rela- 
tion 



A 



fi 2 = m 2 + - (f + (0 2 ) eq .) 



(164) 



Also, from the linearized Vlasov equations we obtain 



U->4 



v p -k)<y/ p = — (v p -k)/? 



2uJ T 



~f 2u p > 



X (l + /oq.(p))/eq.(P) 



O^+En ( S fp' + 6 f-P' + 6 9 P > + fy-p' 



(165) 
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and 



+ 2lu p ) 5g p = — — 



+ E o ( 6 fv' + 6 f~P' + 6 9p' + 5 9- P > 



/eq.(p) 

(166) 



where we have used 



V7 P 

VpWp = — 



V p / eq .(p) = -Vp/3 (1 + / Cq .(p)) /eq.(P) 

Using (11631) . we eliminate 5(f> in (11651) and (11661) and find 



(167) 
(168) 



(u + - V p ■ k)5/ p = ^(V p ■ k)/3 (1 + /eq.(p)) /eq.(p) 



X 1 + 



A 



(cj+ - 2Wp) <5# p = /cq.(p) 1 



C(J^ — k 2 — yU 2 



J 7 



cj 2 — k 2 — yU 2 / 



J 7 



where we have introduced the notation 



P m> 



(169) 
(170) 



(171) 



for the deviation of the fluctuation. Solving (11691) and (j!70j) for 5f p and 5g p 
respectively, and inserting the results into (11711) we obtain the relation 



T = fi(cu+,k)^ 



(172) 



where 



fi(w + ,k) = £ ( 1 + 



A0 2 



x 



d 3 p 1 

cu 2 -k 2 -/i 2 ;7 (2^p^| 

4^ p / cq .(p) | /3(v P -k) 2 
a; 2 , - 4w 2 ^ 2 - (v p ■ kf 



[l + /eq.(p))/eq.(P 



(173) 
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Applying the prescription, 

lim = V — in5(uj - u ) (174) 

e^+o uj — ujq + ie UJ — UJq 



with V implying that the principal part should be taken in integration, we 
find that fl(u + , k) consists of the real and imaginary parts: 

Q(u+,k) = ftx^k) +itt 2 (tu,k) (175) 
The real part is given by 



k) = ~ (l + -P u2 _ X $_f ) *M k) (176) 



with 



4Wp/eq.(p) . /?( V p- k ) 



2 



(l + /eq.(p))/eq.(P) 



uj 2 - Aul u 2 - (v p • k) 5 

(177) 



while the imaginary part is given by 



\2i2 

+— ^tt (<J(o; - cj k ) - 5(cj + cj k )) $i(o; k , fc) 
4cj k 



(178) 



with 



$ 2 (cu, fc) = / [2vr - 2cu p ) - 8{u + 2cu p )) / eq .(p) 

+7T/?V p • k (5(U - V p • k) - 5(UJ + V p ■ k)) (1 + / eq .(p)) /eq.(p)] 

(179) 

where the integration can be carried out analytically, yielding 



•»<". k > - il^r^* - «> + ^^^^ - ^ 

evi-W) — 1 

(180) 
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for u > 0. The values of this function for u < can be found by noting that 
it is an odd function of u. 

The first term in (1180p corresponds to the space-like (scattering) continuum 
with u = ojp+k — uj p ~ v p • k < k. while the second term corresponds to the 
continuum of thermally induced pair creation/annihilation of mesons with en- 
ergy u > 2co>k/2 — 2/i, The function $ 2 (w, k) (and f2 2 (o;, k)) has non-vanishing 
supports in the regions on (u, k) plane as indicated by the shaded areas in Fig. 
2. We plotted in Fig. 3 <&i(a>, k) and $2(^1 as functions of u;/ at = 0.5/x 
for two different values of [ijT. 



2n 



k ' 

Fig. 2. The regions where ft,2(u, k) and &2(u, k) have non- vanishing value is shown 
by shaded areas (two-quasi-particle continua). The solid hyperbola corresponding 
to the meson poles (u = Wk) ° n which ^2(0;, k) has a <5-function singularity in the 
low temperature phase. The shaded area in the time-like region below the curve 
u = 2tJk/2 is only an artifact of the long wavelength approximation. 

The condition 

l=fi(w,k) (181) 

which follows from (11721) determines the dispersion relation of a possible long 
wavelength collective excitation of the system. 

Before examining the solutions of the dispersion relation (I18ip . we note that 
at zero temperature, namely in the absence of quasi-particle excitations, the 
only solution of our coupled kinetic equations is the one which satisfies the 
linearized Klein-Gordon equation (11481) or (11631) with 5f = 5g = and 
/J, = ^Jm? + A0q/2. This gives a simple meson pole u = u k = ^k 2 + jj? 
which appears in the time-like region (u > k). There is no space-like mode 
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- 0.005 ^ ' V ' ■ ■ ■ ' ' ' ' ' ' ' ' 1 

\ 2 4 6,8 

co=k m/k 




Fig. 3. 3>i(cj, k) and ^(^A k) as a function oiuj/k at k/fx = 0.5 and fi/T = 1 (dashed 
line) , jtx/T = 2 (solid line). Both 3>i and $2 vanish at w = fc. $2 also vanishes at 
the pair creation threshold ui = 1[i (in the long wavelength approximation) while 
$1 has a cusp at this point. 



of excitations in the absence of quasi-particle excitations due to the mass gap 
(Au> = 2fi) for the excitations of the vacuum. 

The situation is different in the case of the ordinary Bose-Einstein Conden- 
sate (BEC) which possesses, even at zero temperature, a low-lying acoustic 
mode (Bogoliubov phonon) in the mean- field approximation to repulsive two- 



body contact interaction 411. Il4j. The ordinary BEC contains space-like quasi- 
particle excitations even at absolute zero temperature which corresponds to 
the excitation of one of the zero energy particles forming the condensate. The 
Bogoliubov phonon is a collective excitation of such space-like quasi-particle 
excitation modes modified by the repulsive mean-field interaction. 
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At finite non-zero temperature, the system contains continuum of the space- 
like excitations for all (u, k) satisfying u < k as signified by the non-vanishing 
value of Q,^{uj, k) in addition to the continuum in the time-like region uo > 2fi. 
(See Fig. 2) Any solution of (118 ip in the space-like region therefore is subject 
to the collisionless dissipation known as the Landau damping[39j. 

We plot the function £l 1 (u,\s) in Fig. 4 as a function of u at fixed value 
of k = 0.5 fx. In the low temperature phase, it contains a singularity at the 
position of the mesonic pole at u = \Jk 2 + fi 2 while this singularity disappears 
in the high temperature phase. In the presence of the condensate, a disturbance 
in the quasi-particle distribution may be absorbed into an excitation of the 
condensate which then propagates with the mesonic dispersion relation and is 
converted back again to the quasi-particle excitations. This coupling between 
the condensate and the quasiparticle excitations generates a collective mode 
in the low temperature phase. In the high temperature phase, the quasi- 
particle excitations couple each other only through their direct interaction, 
hence no mixing with single mesonic pole. 

The high temperature behavior of Q\ is the same as that of the function <&i: 
it is negative at u = and increases with u and reaches a positive maximum 
slightly below u = k. In the low temperature phase, Qi(u, k) reverts the sign 
for u below this singularity. One may interpret that effective coupling strength 

*W) = A(l + =r ^r ? )=A(l + ; ^ =7 ) (182) 



changes its sign below the meson pole. Note that here we have used the 
relation A0q = 3fi 2 (111) for the condensate amplitude. We note that Qi(u, k) 
always vanishes at u = k since $i(/c, k) = 0. The cusp at u = 2fi appears at 
the threshold of the two quasiparticle creation. The condition (11811) is fulfilled 
only near the meson pole in the low temperature phase and it gives a shift of 
the mesonic excitation spectrum. We could not find any additional solution 
satisfying ( 118ip . 

We like to note here that without the coupling to the fluctuations to pair cre- 
ation or annihilation, we would not have gotten the first term in the function 
<&i(cj,k) and then Qi(u,k) would have become 1 at u/k > 1 creating a un- 



5 This mode may be compared to the collisionless acoustic mode (the quasi-particle 
sound [48| ) which appears in the superfluid 4 He due to the excitation of the con- 
densate coupled with quasi-particle excitations. 

6 This behavior reminds us of the change of the effective two-body interaction in 
the Bose-Einstein condensates as a function of the external magnetic field due to 
the coupling to the intermediate atomic resonance state, the phenomenon known as 



the Feshbach resonance 14]. 
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damped tachyonic sound mode. Hence the off-diaginal Wigner functions plays 
an important role in making our framework consistent with causality. 
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We introduce the "response function" defined by 



R(u, k) = — Im 



1 — Q(u>, /j) 



fi 2 (w,/i) 



;i83) 



It is plotted in Fig. 5 as a function of uj and compared with the "bare" response 
function given by 



R (uj } k) = -0 2 (o;,/i) 



;i84) 



The response function in the space-like momentum region gives the dynamic 
form factor of the cross section of the scattering of particles coupled to the 
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excitations of the system 4g, |47| |. while its time-like component may give 
the rate of pair annihilation of the quasi-particles. We note that in the low 
temperature phase the response function changes its sign below the meson 
pole due to the sign change of the effective coupling (11821) . We observe some 



enhancement (depletion) of the strength in the time-like pair annihilation near 
threshold and small depletion (enhancement ) in the space-like region in the 
low (high) temperature phase. 
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Fig. 5. Plot of the response function R(u,k) (|184|) (solid lines) at kj\i = 0.5. For 
comparison the bare response function RQ(ui,k) = —^{uj^k) is also shown (dashed 
lines). Note that in the low temperature phase R changes sign in the space-like 
region and at the meson pole (to ~ tu^) due to the sign change of the effective 
coupling strength (1182j) , 

As we have noted, in the low temperature phase with non-vanishing con- 
densate amplitude (fio, the mesonic excitations couple with the quasi-particle 
excitations and this give the additional shift of the meson mass from the one 
obtained by the gap equation. The shift of the meson pole from the bare 
spectrum u = uj^ = \/k 2 + n 2 may be computed by the condition 



2 ; 2 2 

uj — k — fi 



k 2 - fi 2 )Q(uj, k) = 



(185) 
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which in the long wavelength limit ( k = ) gives a solution at u 2 = fi' 2 . The 
shift of the meson mass Afi 2 = fi' 2 — fi 2 is determined by 



where we have used A0 2 , — 3/i 2 . In the case of |A/i 2 | << /i 2 this may be solved 
approximately 



The equation (I186p may be solved by iteration starting from the first approx- 
imation (I187P by inserting it for A/i 2 in the integrand on the right hand side. 
We plot in Fig. 6 the meson mass shift determined by this method. The me- 
son mass becomes smaller as T increases and eventually becomes zero at the 
temperature which satisfies the condition 



1 = A /(l^4 /o, '( p) (188) 

The vanishing of the effective mass of mesonic excitation may be identified 
as the onset of instability in the metastable low temperature phase by small 
fluctuation (spinodal decomposition). The same instability occurs when one 
approaches to T = T c from high temperature phase. 



We note that the temperature T sp at 
not coincide with the " backbending" 
The latter temperature is determined 

1 = A /(^2^ /e - (p) 



the spinodal point given by (11881) does 
temperature T\ where dfi/dT diverges, 
by the condition 

(189) 



which can be obtained from the gap equation flllOD . We found that T\ is 
slightly above T sp and the unstable region appears associated with the lower 
solutions fi of the gap equation. These unstable solutions thus appear only in 
the region which is not easily accessible and may well be considered as another 
pathology of the mean field approximation. 



7 These two temperatures may coincide with T c in the case of the second order tran- 
sition of the Landau-type. Such behavior has been obtained in 36f] by a temperature- 
dependent loop expansion method formulated in [351 ] . 
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Fig. 6. Mass shift of the mesonic excitation: solid line (//) is the mass of the exci- 
tation determined by (|186p and the dashed line is the solution of the gap equa- 
tion evaluated at A =10. The slopes of two curves diverge at the temperature 
t\ = T\/T c = 1.0714 while fi' vanishes at t sp = T sp /T c = 1.0702. 

7 Concluding remarks 



In this paper, we have developed a kinetic theory for a system of interact- 
ing quantum fields in the mean field approximation taking into account the 
existence of quasi-particle excitations. We have obtained a set of coupled equa- 
tions of motion, one for the meson condensate in the form of non-linear Klein- 
Gordon equation which contains extra non-linearlity due to the particle exci- 
tations. The equations of motion for quasi-particle excitations are described in 
terms of the Wigner functions, which reduce to a semiclassical Vlasov equation 
for one-particle distribution function with a modification due to the coherent 
pair-creation and pair-annihilation expressed by the off-diagonal components 
in the Wigner function. These off-diagonal components may be eliminated by 
a suitable Bogoliubov transformation of the particle creation and annihila- 
tion operators for a uniform, time-independent system. However, they remain 
non-vanishing in general non-uniform, time-dependent systems. 



We have shown that in equilibrium these equations are reduced to a gap 
equation in the Hartree approximation. This implies that our kinetic equations 
are natural extension of the Hartree approximation to the non-equilibrium 
situation. It is well known, however, that in this approximation the phase 
transition becomes first order. This is a generic feature of the mean field 
approximation 27, 37, 34| which persists also for multi-components scalar 



fields with continuous O(N) symmetry for arbitrary finite integer N. 

In the present work, we studied also the excitations spectrum in the system 
near equilibrium. In the high temperature phase there is no collective ex- 
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citation mode in the system besides the two continua of the quasi-particle 
excitations in the entire space-like energy-momentum region and the time-like 
region above the pair creation threshold. In the low temperature phase we 
found that coupling of the meson pole to the quasi-particle continua give rise 
to the shift of the meson mass which becomes zero at the edge of the spinodal 
instability line of the first order transition. The spinodal point appears deep 
inside the meta-stable region of the first order transition which may well be an 
artifact of the mean field approximation. We found that the coupling to the 
off-diagonal components of the Wigner functions should be properly included 
to avoid the appearance of the undamped tachyonic mode. 



In this work we used a single component real scalar field model which possesses 
only discrete symmetry. It is straightforward to extend the present analysis 
to models with continuous symmetry such as the sigma model with O(N) 
symmetry. Basic features of the present analysis is preserved in such extension. 
It is known however that the Goldstone theorem is apparently violated in 
this approximation for a system with continuous symmetry. We will show in 
the forthcoming paper 22j, that the missing Nambu- Goldstone mode may be 
retrieved in the collective excitations of the system. 



Our coupled kinetic equations may be solved for an arbitrary initial conditions. 
We plan to study the freezeout dynamics with these equations with more 
realistic interactions. It would be interesting to see in particular how much flow 
is generated by the acceleration by the mean field as the vacuum condensate 
is restored. 
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A Proof of the equivalence of the Gaussian density matrix average 
and the mean field approximation. 



Here we present some details of the computation of the equation of motion 
of the Wigner functions and show that with the Gaussian Ansatz for the 
density matrix the result is equivalent to what we obtain from the mean-field 
Hamiltonian 



We like to compute the time-derivative of the operator product a pi a p2 which 
appears in the definition of F(p, k, t) with pi = p + k/2 and p 2 = p — k/2. 



■ d ( t 



P2 



ia pi a P2 + ia pi a P2 
[a pi ,#]a P2 + a p Ja P2 ,#] 



(A.1) 



The commutators [a pi , H] is decomposed into a sum of the commutators with 
Hi ( i — 1, • • • , 4) among which the commutators with Hi and H% do not sur- 
vive the average with the Gaussian density matrix since they only contain the 
odd power of the field operators. We only need to compute the commutators 
with H 2 and H^. 

To compute the commutator with H 2 it is convenient to rewrite H 2 as 

H 2 = £ u p (a p a p + \) + ^J drtfir, t)0 2 (r, t) (A.2) 

We then find 



[a p (t),H 2 ]=uj p a p (t) + - I dr#(r,t)[a p (t),0 2 (r,t)] 



A 



-ipr 



■u p a p {t) + dr<j) 2 c (r,t)^=(j)(r,t) 
aUt),H 2 }=-u p aUt) + ^ I dv<Pl(r,t)[aUt)^ 2 (r } t)] 



= -uj p a p {t) -~j c/r0 2 (r, t )^^= ( 
where we have used the following formulae: 



(A.3) 



(A.4) 



"p' 



■ n< 



t.n-1. 



r,t)k,fet)] 



-n> 



Ln-l, 



T,t) 



ipr 



(A.5) 
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"pi 



r,t) =n0"- 1 (r,t)[a p ,</.(r,t)] =n0^ i (r,t) 



V 7 ^ 



With these results we obtain 



W Pl a P 2, H 2 ] = -(w pi - cu P2 )a pl a P2 



~2 7 ^( r ,t)0(r,t)-^==a P2 



ipi-r 



+ - 



| dr<f> 2 c (r,t)4>(r,t) 



The Gaussian average of this equation gives 



([a Pl a P2 , // 2 ]) = -(w pi - w P2 )(a Pi a P2 ) 



A 



3 ipi-r 



~2 y dr ^c( r ^)^=(0( r >o° P 2> 



+- 



A 



/ g -iP2-r 
dr^( r ,*)— = (^(r,*)a 



Pl/ 



Next we compute the commutators with H 4 . We first compute 



A r ~ e ipr 
[a p (t),H 4 } = -j dr^(r,t)^=a p (t) 



,ipr 



f p (t),H 4 ]=-± J drft(r,t)-^=a p (t) 



We then use these results to obtain 



A f 



o«Pi-r 



*T2 



Pl 



A 



-iP2-r 



+- / dr^(r,i)- 



pi 



Taking the average with the Gaussian density matrix, we find 



([a Pi a P2 ,H 4 }) = -- J dr((j) 2 (r,t)) —= {<j)(r,t)a P2 ) 



ipi-r 



+ - 



A 



/ p -ip2-r 
dr(^(r,t))— =<0M)at 
V ZUJ P2 
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We observe the similarity between the last two terms of the commutator with 
H 2 and the commutator with if 4 after Gaussian average. These terms can be 
combined by introducing the self-energy function, 

II(r,f) = #(r,f) + (0 2 (r,t)> = £n q (t)e^ r (A.13) 

q 



Integration over the space coordinate r give a delta function <5(q + p + pi) 
and 5(q + p — P2). Performing the integral over p yields 



( [4l G P2 > H ] ) = - ( W Pl - U P2 ) (Opi fl P2 ) 



A E / (( a -pi-q a P2) + ( a pi+q a p2)) 

4 q v^pi+q^ps 

+ T E 7= (( Q P2-qO + (a- P2 +qO) 

4 q V W Pl W P2-q 



(A.14) 



This result coincides with the average of the commutator with the mean field 
Hamiltonian H m { defined by (l68j) . 

([a pl a P2 ,iJ]) = ([o^Opj,,^]) (A.15) 



Commutators of the four bilinear operator products with the mean-field Hamil- 
tonian are listed below: 



[at. ctr.^ , H m A 



Pl P2' 



(w pi - ^ P2 )a Pi a P2 



4e 



q V^Pl+q^P 2 



1 x 



2 q V^Pi^Pa-q 



_ Pl _ q a P2 + a P1+q a P2 ) 
(a pi a P2 _ q + fl Pl fl- P2 +q) 



[a pi a P2 , H m {] — (a? pi + ^ P2 )a Pl a P2 



n r 



q \/^Pl+q a; P2 



(a Pl _ q + a_ Pl ^q)a, 



P2 



1 . - 



2 q ^Wp^pa-q 



a pi( a P2-q 4~ a -p2+q) 



a Pi a p2 ,H, 



mf 



-(u pi +^ P2 )a Pi a P2 



-~E 

2 q v /a; Pl+q ci;p 2 



v a -pi-q a P2 "I" a pi+q a p 2 



(A.16) 



(A.17) 
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^ — — ( a P2-q fl pi + a -p2+q a pi, 



2 q V^Pi^Pa-q 



(A.18) 



a Pl a P2 , H mi = (u> pi - cu P2 )a Pl al 



P2 



+ z2 



"pi-q 



l -pi+q) a P2 



2 q V^W^W-q 



a p! ( a P2-q + a -p2+q) 



(A.19) 
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